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ABSTRACT 

In this paper, we introduce definition for fractional 
quadruple Laplace transform of order a,0 < a < 1, for 
fractional differentiable functions. Some main 
properties and inversion theorem of fractional 
quadruple Laplace transform are established. Further, 
the connection between fractional quadruple Laplace 
transform and fractional Sumudu transform are 
presented. 

KEYWORD: quadruple Laplace transform, Sumudu 
transform, fractional Difference. 

INTRODUCTION 

There are different integral transforms in mathematics 
which are used in astronomy, physics and also in 
engineering. The integral transforms were vastly 
applied to obtain the solution of differential equations; 
therefore there are different kinds of integral 
transforms like Mellin, Laplace, and Fourier and so 
on. Partial differential equations are considered one of 
the most significant topics in mathematics and others. 
There are no general methods for solve these 
equations. However, integral transform method is one 
of the most familiar methods in order to get the 
solution of partial differential equations [1, 2]. In [3, 
9] quadruple Laplace transform and Sumudu 
transforms were used to solve wave and Poisson 
equations. Moreover the relation between them and 
their applications to differential equations have been 
determined and studied by [5, 6]. In this study we 
focus on quadruple integral determined and studied by 
[5, 6]. In this study we focus on quadruple integral 
transforms. First of all, we start to recall the definition 
of quadruple Laplace transform as follows. 


Definition: Let / be a continuous function of four 
variables; then the quadruple Laplace transform of 
f(w, x, y, z ) is defined by 


-‘wxyz 


f(w,x,y,z ) = F(p,q,r,s ) 


( 1 ) 
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Where w,x,y,z> 0 and p,q,r,sare Laplace 
variables, and 

f(w,x, y, z ) 
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is the inverse Laplace transform. 

Fractional Derivative via Fractional Difference 
Definition: let git) be a continuous function, but not 
necessarily differentiable function, then the forward 
operator FW( h) is defined as follows FW(h)g(t) — 
g(t + h ), Where h>0 denote a constant 
discretization span. 

Moreover, the fractional difference of g(t) is known 
as 

A a g(t) — (FW — h) a g(t ) 

oo 

= ^ (-l) m Q git + (a - m)h ) 

m =0 
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Where 0 < a < 1, 

And the a-derivative ofg'(t) is known as 
, . A a g(t) 

„(a)( t ) _ -- 

y w hlo h a 
See the details in [9, 10]. 


Fractional quadruple Sumudu Transform 

Recently, in [13] the author defined quadruple 
sumudu transform of the function depended on two 
variables. Analogously, fractional quadruple sumudu 
transform was defined and some properties were 
given as the following 

Modified Fractional Riemann-Liouville Derivative 

The author in [10] proposed an alternative definition Definition: [14] The fractional quadruple sumudu 
of the Riemann-Liouville derivative 


Definition: let g(t ) be a continuous function, but not 
necessarily differentiable function, then 
Let us presume than git) — K , where K is a 
constant, thus a-derivative of the function g(t) 

isD?* = [ '« S0 ' 

1 l0, otherwise. 

On the other hand, when g(t) =£ K hence 
git) = g(0) + (git) - g(0)), 

and fractional derivative of the function git) will 
become known as 

g (a \t) = D?gi0) + D?(git) - p(0)), 
at any negative a, (a < 0) one has 
D?(git) -giO)) 

t 

= r( _ a ) Jit- gig)dg, a < 0, 

o 

While for positive a, we will put 

O“GK0 ~9(0)) = D‘g(t) = D t (g(“-0). 

When m < a < m + 1, we place 

g (a) it) = ig ia - m) it))^ m \ 

m < a < m + 1 ,m > 1. 

Integral with Respect to (dt)“ 

The next lemma show the solution of fractional 
differential equation 


dy — git)(dt) a ,t > 0,y(0) = 0 
By integration with respect to ( dt) a 


( 2 ) 


Lemma: If git) is a continuous function, so the 

solution of (2) is defined as the following 
i 

yit) = J giv)idrj) a , y(0) = 0 


L 

— a j it — g) a ~ 1 gig)dg, 0 < a < 1 


( 3 ) 


For more results and varies views on fractional 
calculus, see for example [15, 16, 17, 18, 19, 20, 21] 


transform of function f(x, t) is known as 

Sa{fix,y,z,t)j = Gaip,q,r, s), 

n GO n GO n CO n CO 

= J J J J E ai~ix + y + Z 

+ t) a )f(px, qy,rz, st)idx) a idy) a idz) a idt) a 


Wherep, q,r,s G C, x,y,z,t > 0 


and 


Z co 

i 

m= qJ 


X 


t = 0 r(am + 1) 
is the Mittag-Leffler function. 

Some Properties of Fractional Quadruple Sumudu 
Transform 

We recall some properties of Fractional Quadruple 
Sumudu Transform 

fiax, by, cz, dt)} — G^iap, bq, cr, ds), 

S%{ fix — a,y — b,z — c,t — d)} 

— E a (—(a + b + c + d) a )Gaiap, bq, cr, ds), 
Sa{d?fiax, by, cz, dt)} 

Ga iv> q,r,s) - Til + a)/(0, x, y, z) 


u u 


(4) 


where df is denoted to fractional partial derivative of 
order a (see [10]). 

Main Results 

The main results in this work are present in the 
following sections 

Quadruple Laplace Transform of Fractional order 
Definition 6: If f(x,y,z,t) is a function where 
x,y,z,t > 0, then Quadruple Laplace Transform of 
Fractional order of f(x, y, z, t) is defined as 
Laifix, y> z,t)} = F^ip, q, r, s) (5) 

/»GO /»GO /'GO n GO 

= Jill E ai~ipx + qy + rz 
+st) a )fix,y,z, t)idx) a idy) a idz) a idt) a 

Where p,q,r,s E C and E a ix) is the Mittag-Leffler 
function. 
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1 . /p q r S' 


Corollary 10.1 By using the Mittag-Leffler property 
then we can rewrite the formula (5) as the following: 
E 4 a {f(x, y, z, t)} = Fa(p, q, r, s') 

p OO p 00 p GO p GO 

= J J J J E a(-(P x ) a ) E a(-qy a )E a (-(rz) a ) 

E a (-( St) a )f{x, y, z, t) ( dx) a (dy) a (dz) a (dt) a (6) 

Remark 8: In particular case, fractional quadruple 
Laplace transform (5) turns to quadruple Laplace 
transform (1) when a = 1. 


—_ ± 

a a b a c a d a a \a’b’c’dJ 


3. Shifting property 

Let L 4 a {f (x, y, z,t)} = F 4 ( p, q, r, s) then 

Ea[E a (—(px + qy + rz + st) a )f(x,y, z, t)} 

= F^ (p + a,q + b,r + c,s + d) 

Proof: 

L, 4 a {E a (—(cix + by + cz + dt) a )f(x, y, z, t)} 


pCfJ p ou p ou p ou 

= E a (-(ax + by 

j 0 Jo Jo Jo 


+ cz 


Some properties of fractional quadruple Laplace 
transform 

In this section, various properties of fractional 

quadruple Laplace transform are discussed and proved By using the equality 
such as linearity property, change of scale property 
and so on. 


+ dt) a )E a (—(px + qy + rz 
+ st) a ) f(x,y,z,t)(dx) a (dy) a (dz) a (dt) a 


[E a (\(x + y + z + t)“) 
=E a (Ax a ) E a (Ay a )E a (Az a )E a (At a ) 


1. Linearity property 

Let (x, y, z, t) and f 2 (x, y, z, t) be functions of the 
variables x and t, then 
Eaia-Mx.y.z, t) + a 2 f 2 (x,y,z, t)} 

=a i E 4 a {f i(x, y, z, t)} + cl 2 L 4 a {f 2 (x, y, z, t)}(7) 


Where a l and a 2 are constants. 

Proof: We can simply get the proof by applying the Hence 
definition. 


Which implies that 

La{E a (—(cix + by + cz + dt) a )f(x,y,z, t) 

p 00 p CO p GO p GO 

= J J J J E a (~((a + p)x + (b + q)y + (c 

+ r)(d + s)t) a ) f(x,y,z, t)(dx) a (dy) a (dz) a (dt) a 


hence 


2. Changing of scale property 
If l 4 a {f (x, y, z, t)} = F 4 (p, q, r, s) 

L 4 a {f(ax, by, cz, dt )} 

= _i_ 

a. a b a c a d a a b c d' 


Whenevera and b are constants. 

Proof: 

{/(ax, by, cz, dt)} 

p GO p GO p 00 p GO 

= 1 I I I E a (-(px + qy + rz 
+ st) a )f(ax, by, cz, dt) (dx) a (dy) a (dz) a (dt) a ( 8) 

We set j — ax , k — by, l — cz and m — dt in the 
above equality (8), therefore we obtain 
L 4 a {f(ax, by, cz, dt)} 

1 rr r r E (_(Pi + ^ + r _i 

a a b a c a d a J 0 J 0 J 0 J 0 “\ la b c 
+ —J J f(J, k, l,m)(dx) a (dy) a (dz) a ( < dt) a 


L a {E a (-(px + qy + rz + st) a )f(x,y,z, t)} 
=F 4 {p + a,q + b,r + c,s + d) 

4. Multiplication by x a t a 
Laif(x, y, z,t)} = F 4 {p, q, r, s) 


p GU p OU p GU p OU 

= J J J J E a (-(px) a )E a (-(qy) a )E a (-(rz) a ) 

E a (-(st) a )f(x, y, Z, t)(dx) a (dy) a (dz) a (dt) a 
then 

d 2a 

[x a y a z a t a f(x,y,z,t)} — 


E 4 a {f(x,y, z, t)} 


dp a dq a dr a ds a 


Proof: 


L 4 a {x a y a z a t a f(x, y, z, t )} 


p GO p GO p GO p GO 

= J J J J x a E a (—(px) a )E a (—(qy) a )E a (—(rz) a ) 

E a (—(st) a )f(x,y,z,t)(dx) a (dy) a (dz) a (dt) a 
By using the fact D? (E„(— s“x“)) = 
x a £' a (—s a x a ), then 


L 4 a {x a y a z a t a f(x, y, z, t )} 
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r cc r co r co r co d a d a Proof: We obtain the definition of fractional 

J J J J Qpu ^a\ (P x ) ' Qqa^ay (MV) ) quadruple Laplace transform and fractional quadruple 
Q a g a convolution above, then we obtain 

— E a (-(rz) a — E a (-(st) a )/(x, y, z, t)(dx) a (dy ) a g )( x , y, z, t)} 

y' oq oq oq ^ r°° r 00 r 00 r 00 

= [ /„ /„ l dp a dq a dr a ds a j 0 J„ J„ J„ £ “( _( P* ) “) £ »(“ (d >' ) “) £ »(“ (rz) “: ) 

Ea(~(. rz ) a )E<x(—(st) a )f(x,y, z, t)(dx) a (dy) a (dz) a (dt^^“^p a 2^^ ****« p)(dx) a (dy) a (dz) a (dt) a 

= [ \ f [ E a (—(px') a> )E a (—(ciy) a> )E a (—(j'z) a> )[dt) a 
j 0 ‘'O ‘'O ‘'O 


Remark: All results above are suitable for quadruple 
Laplace transform when a = 1. 


£ a (-OO a ) 


f* f z 

f(x — rj,y — 9,z — y,t 

J o Jo Jo Jo 

(dx) a 


- d)#07,0, y, 5)(dq) a (d0) a (dy) a (d5) a 
(dy) a (dz) a (dt) a 


Theorem: If the fractional quadruple Laplace 
transform of the function /(x, y, z, t) is 

LUnx,y,z,t)} = F^p,q,r,s), and fractional i et ; = x - rj, k = y - 9, l = z - y, m = t - 8 and 
quadruple sumudu transform of the function 
/(x,y,z,t) is 5^{/(x,y,z,t)} = G«(p,q,r,s),then 

„ 1 /I 1 1 1 

Ga(.P> q,r,s) - a E a (->->->- 


pu.qu. r u. S L 


p q r s 


Proof: 

Ga(p> q>r,s) 

n GO /'CO /• CO /• 00 

-J J J J Ea(-x a )E a (-y a )E a (-z a )E a (-t a ) 

fipx, qy, rz,5t)(dx) a (dy) a (dz) a (dt)“ 

By using change of variables _/ -> px, k -» qy, Z -» rz 
and m d 


taking the limit from 0 to oo, it gives 

p CO /»CO p CO p CO 

= J J J J Ea(-p a 0 + ri ) a ) 

E a (-q«(k + 9) a )E a (—r a (l + Y) a )E a (—s a (m 

-«//// /(/, k, l, rri)g(rj, 9, y, 5) (dq) a 
(d9) a (dy) a (d8) a (dj) a (dk) a (dl) a (dm) a 


n OO p- CO /•OO /'OO 

= J J J J E a {-p a nE a {-q a k^E a {-r a l a 

E a (.—s a m a )f(J, k, l, m)(dj) a (dk) a (dl) a (dni) a 
x / rv ,fW 


) 


' tvM" f f f '• (-©") '• (- (;)') 'L(f.( S | 

/ K(-s a $ng(v,e,y,5)<id V rm a (dYy(d6y 

E a ( — y —J \fQ,k,l,mydjy(jdky(jdiy(jdm) a 

A 1 ,/l 1 1 h 

G a (.P>q> r > s ) — a a r a s ct ^ a \p'q’r’~s) 


Laifix, y, z, t)}L’{/(x, y, z, t)} 

Inversion formula of Quadruple Fractional 
Laplace’s transform 

The Convolution Theorem of the Fractional Firs Hy. we wiU set U P definition of fractional delta 
Double Laplace Transform func(ion of tw0 variables as follows 

Theorem: The double convolution of order a of 

functions/(x,y,z, t) and p(x,y,z,t) can be defined Definition: Two variables delta function d a (x 
as the expression 
(/(x, y, z , t) * * * * a flf (x, y, z, t)) 

p* py pz pt 

= J J J I fix — q,y — 9,z — y,t 

- 8)gip,9,y,°8) (dp) a (d0) a (dy) a (dd) a , 
therefore one has the equality 
£«{(/ ****« g)ix,y,z, t)} 

= £«{/(*, y, z, t)K4{/(x, y, z, t)} 


a,y — b,z — c,t — d) of fractional order a, 0 < a <1, 
can be defined as next formula 

I I gix,y,z,t)8 a ix - a,y - b,z - c,t 
Jr Jr Jr J 

— d)(dx) a (dy) a (dz) a (dt) a 
= a 4 g(a,b, c, d)(9) 
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In special case, we have 

f f f [ dl x , y> z, t) 8 a (x, y, z, t)(dx) a (dy) a (dz) a (dt) a 
Jr Jr Jr J 


SIS S Ea ^ ix<X(<h - vY)Ea(iy a( <u - q) a ) 


— 


a* g (0,0,0,0) 


E a (iz a (v — r) a E a (it a (w — s) a 
(dp) a (dq) a (dr) a (ds) a 


Example: we can obtain fractional quadruple Laplace f f f f,, 

transform of function S a U - a,y- b,z- c.t-d ) = J J J J (*0»(dy)«(dz)«(dt)“ 

as follows r r 

L\{8 a (x-a,y-b,z-c,t-d )} j j j J E a (i(-xj) a )E a (i(-yk) a ) 

= r r r [\c-cpx+ qy +r Z+s m 

8 a (x-a,y -b,z-c,t- d)(d*)«(dy)« =f III E alx,y,z,t)(dxy(dyy(dzy(dty 

(dz) a (dt) a 

= a 4 E a (—(pa + qy + rz + st)“)(10) 

In particular, we have La{8 a (x,y,z,t)} — a 2 


= a 

Note that one has as well 

„4 


Ty* fill Ea ^~ hx ^ aEa ^~ uy ^ aEa ^~ vz ^ a 


Relationship between Two Variables Delta 
Function of Order a and Mittag-Leffler Function Wa ^ R R R 

The relationship between E a (x + y + z + t) a and E a (i(-wt) a (dh) a (du) a (dv) a (dw) a = 8 a (x,y,z, t) 
8 a ( x > y< z < t) is clarified by the following theorem 

Inversion Theorem of Quadruple Fractional 
Theorem: The following formula holds Laplace Transform 

Theorem: Here we recall the fractional quadruple 
Laplace transform (5) for convenience 
La{f(x,y,z, t)} = Fa(p,q,r, s)(12) 

n 00 n 00 z' 00 n OO 

= E a (-(px +qy+ rz 

J 0 Jo Jo Jo 

+ st) a )f(x, y, z, t)dx a dy a dz a dt a 

And its inverse formula define as 
f(x, y, z, t) 


fill Ea ^~ hx ^ aEa ^~ uy ^ aEa ^~ vz ^ a 

E E E E 

E a (i(-wt) a (dh) a (du) a (dv) a (dw) a 
= 8 a (x,y, z, t)(ll) 

Wher eM a satisfy the equivalence E a (i((M a ) a ) =1, 
and it is called period ofthe Mittag-Leffler function. 
Proof: We test that (11) agreement with 

a 2 


= j J J J E a (i(hx) a E a (i(uy) a E a (i(vz) a 

E E E E 

E a (i(wt) a (dx) a (dy) a (dz) a (dt) a 

We replace 8 a (x,y,z,t) in above equality by (11) to 
get 

“ 2 = I If /(*)Vy)“(*0“(dt)“(£j;; 

E E E E 

ISIS E a (i(hx) a )E a (i(uy) a )E a (i(vz) a ) 

E E E E 

E a (i(wt) a )E a (i(-px) a E a (i(-qy) a E a (i(-rz) a 

E a (i(—st) a (dp) a (dq) a (dr) a (ds) a 

= j j I f(dx) a (dyy(dzy(dt)‘‘-?E s 


■j /'CO /• 00 /• 00 /'CO 

= 777^ E a((px + qy + rz 

V IV ‘cc) Jo Jo Jo Jo 

+ st) a )Fa(p, q, r, s)(dp) a (dq) a ( dr) a (ds) a 


(13) 

Proof: 

Substituting (12) into (13) and using the formula (11), 
(9) respectively, we obtain in turn 


+ loo +loo +loo +100 


fix, 


■ y ’ z ’ t)= dyl I I I EApx) “ 


— ico —ico —ico —ico 


E a (qy) a E a (rz) a E a (st) a (dp) a (dq) a (dr) a (ds) a 

z' CO s' 00 /> CO z' GO 

Jill E a (-(pj+ qk + rl + sm) a 
f(j, k, l,m)(dj) a (dk) a (dl) a (dm) a 
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OO 00 00 00 

= (jj III f(fi,y,^S)(dfir(dyr 

0 0 0 0 

+ioo +ioo +io o +£oo 


o dxpym a Jill E a (P a (x-j) a ) 
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